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We discuss how future cluster surveys can constrain cosmological parameters with particular 
reference to the properties of the dark energy component responsible for the observed acceleration 
of the universe by probing the evolution of the surface density of clusters as a function of redshift. 
We explain how the abundance of clusters selected using their Sunyaev-Zel'dovich effect can be 
computed as a function of the observed flux and redshift taking into account observational effects 
due to a finite beam-size. By constructing an idealized set of simulated observations for a fiducial 
model, we forecast the likely constraints that might be possible for a variety of proposed surveys 
which are assumed to be flux limited. We find that Sunyaev-Zel'dovich cluster surveys can provide 
vital complementary information to those expected from surveys for supernovae. We analyse the 
' impact of statistical and systematic uncertainties and find that they only slightly limit our ability 

to constrain the equation of state of the dark energy component. 
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I. INTRODUCTION 



Sunyaev and Zel'dovich [1,2] (SZ) first noted that cosmic microwave background (CMB) photons would be re- 
scattered as they pass through the hot intergalactic medium in clusters of galaxies due to inverse Compton scattering. 
Since this process preserves the overall number of photons, one observes a decrement in CMB temperature in the 
iy~j ■ Rayleigh-Jeans part of the spectrum, and a corresponding increase at high frequencies above the null frequency of 
217 GHz. For a cluster of given mass, the brightness temperature depends only on the integrated pressure of the gas 
in the cluster along the line of sight and not on its distance since the red-shifting of the photons is exactly balanced 
by their higher density at the time of scattering. This not only allows the clusters to be detected at much higher 
redshifts than using standard optical and X-ray measurements, it also makes any detection much less susceptible to 
fS i 1 the internal dynamics of the cluster, providing a reliable method for detecting clusters in a blank field survey. 

To date the SZ effect has been mapped in a small number of clusters with targeted observations using either 
■ a single dish [3-5], or an interferometer [6-15]. In the very near future blank field surveys will be performed on 
dedicated instruments with arcminute resolution which should yield flux limited samples of clusters at a variety of 
frequencies and resolutions over large areas (see, for example, [16-28] and references therein). Moreover, instruments 
originally designed to make small scale CMB anisotropy measurements such as the Very Small Array (VSA) [29], 
Cosmic Background Imager (CBI) [30] and those on board the Planck Surveyor [31,32] will also inadvertently provide 
information, albeit less efficiently due to their large beam sizes [33], though the Planck Surveyor can balance this 
disadvantage with its full sky coverage and the range of frequencies available. 

The limiting mass of a survey as a function of redshift can be computed from the limiting flux assuming the cluster 
is virialized and, hence, that there exists a relationship between the mass of the cluster and the temperature of the 
gas within it. Therefore, the survey yield is a calculable function of the cosmological parameters and the process 
of halo formation. If we believe the physics of halo formation is understood, which in the case of cluster size halos 
should just involve gravitational physics, then the number of clusters as a function of redshift can be used to constrain 
cosmological parameters [34-45]. More precisely the number of objects will depend on the comoving volume element, 
the angular diameter distance and the rate of growth of perturbations, all of which depend sensitively on the late-time 
evolution of the cosmological scale factor a(t). 

The study of the late-time evolution of the universe has become an important part of many observational pro- 
grams ever since observations of type la supernovae (SNe) suggested that the expansion of the universe might be 
accelerating [46-49]. It was found that a combined initial sample of 42 SNe with (z) ~ 0.6 were much dimmer at 
high redshift than required by a matter dominated Einstcin-de Sitter universe, with the simplest remedy being the 
inclusion of Einstein's 'greatest blunder', the cosmological constant A which gives rise to acceleration by virtue of the 
negative pressure associated with vacuum energy. Subsequent observations of more SNe have only strengthened this 
conclusion. These observations probe the magnitude-redshift relation by assuming the SNe are standard candles. 

However, this is notoriously difficult to incorporate into a realistic particle physics theory since estimates of the 
vacuum energy, while very rarely zero, are substantially larger than the observed value. This lead to the proposal of 
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a dynamical dark energy component with negative pressure, often called 'Quintessence', due to a scalar field [50-62]. 
Although this idea is not without its own problems, it provides an interesting theoretical framework to parameterize 
and test the idea of dark energy. 

If one takes Quintessence seriously one has to accept the fact that the equation of state parameter wq = pq/pq 
is not only negative, but that it can evolve with time. There are many specific dark energy models, none of which 
are theoretically compelling and, therefore, in order to test them one requires some kind of parameterization of this 
evolution. There has been some debate in the literature about what is the best way to do this: suffice to say that it 
is difficult to imagine one which accurately represents every model and it is often best to tailor the parameterization 
to the type of observations under consideration [63-71]. Since the clusters under discussion here will have redshift of 
less than 1.5, a Taylor expansion of the form wq = wo + W\z + .. , where Wq and w\ are constants is a sensible form 
to consider. 

An important feature of SZ surveys, and similar observations that probe number counts, is that their dependence 
on the expansion rate of the universe is very different to that of SNe observations, that is, they are complimentary 1 
and, hence, can be used to break parameter degeneracies. A degeneracy between Wq and w\ is present for both SZ 
and SNe surveys. However, the shape of the likelihood contours in the wq-w\ plane is very different [45], allowing one 
to use the two together to make more substantial statements as to the nature of the dark energy. 

In this paper we will discuss using the evolution of the surface density as a a function of redshift. A particular 
important ingredient for using number counts of clusters to constrain cosmology is the ability to measure their redshift. 
Ideally spectroscopic redshifts (Az ~ 0.01) would be available for each cluster. However, it may only be possible to get 
photometric redshift (Az s=a 0.1). Moreover, at present it is difficult to estimate redshifts for objects with z <~ 1. The 
Sloan Digital Sky Survey (SDSS) [74,75] and the Visible and Infrared Survey Telescope for Astronomy (VISTA) [76] 
will measure the redshift of many hundreds of thousands of galaxies and will provide this information at least to some 
degree. Except where explicitly stated we will use Az = 0.1 and z max = 1-5, optimistically assuming techniques to 
find redshifts for very distant objects can be developed on a similar timescale to instruments under discussion here. 

In section II we will discuss the principal cosmology dependence of the surface density of clusters on cosmology 
before going on to discuss the relation between the limiting mass and the survey set up in section III. We will then 
perform a mock likelihood analysis for six proposed survey setups in section IV and discuss statistical and systematic 
uncertainties in our assumptions in section V. Note that throughout this paper algebraic relations will be expressed 
in natural units where c = h = = 1 • 



II. SURFACE DENSITY OF CLUSTERS 



There exists a plethora of dark energy models none of which are particularly compelling. We will only discuss 
models which are minimally coupled to gravity. Since, a priori, there is no theoretical preference for any particular 
model we take a phenomenological approach and characterize the dark energy component by a linear evolving equation 
of state [66-68] 

wq = w + w\z . (1) 

This is sufficient since we are only interested in the low redshift behaviour of the dark energy component and further- 
more (1) approximates a wide range of dark energy models adequately at low redshifts [66-69]. 

In order to predict the surface density of clusters of a mass limited Sunyaev-Zel'dovich survey we need to compute 
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where dV/(dzd£i) = [r(z)] 2 /H is the comoving volume in a flat universe, with r(z) = J* H~ 1 (z')dz' the coordinate 
distance, Afl is the angular sky coverage of the survey and dn/ dM is the comoving number density of objects with 
mass M and redshift z, sometimes called the mass function. Mi im is the limiting mass, which will depend in general 
on the parameters of the survey. In this section, we will assume that Mu m is constant and does not depend on the 



1 Note that measurements of the angular diameter distance relation at low redshifts using, for example, gravitational lensing 
statistics [72,73] have exactly the same degeneracy as those using the magnitude-redshift relation. These can be used as a 
direct test of the veracity of SNe surveys. 
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cosmological parameters. In the next section, however, we will drop this assumption and model Mn m in a more 
realistic way for our mock likelihood analysis. 

We will use the comoving number density calibrated using a series of N-body simulations performed by the VIRGO 
consortium [77], with 

*L (2 ,M)=-0. 2 2^M^exp{- l 0.73-lo g[ D( 2 ). M ]|- } , (3) 

where the mass, M of the object is defined to be that inside a spherical over-density of A = 200 times the critical 
background density, that is, M 2 oo = 47ri?2o 200/9 cr it(;z)/3. D(z) is the growth factor normalized to have D(0) = 1. 
We obtain the growth factor by solving numerically the perturbation equation for the matter fluctuations 

C + fa" 1 [1 - Ma) (1 - O m (a))] 5' m - ^a- 2 n m (a)5 m = , (4) 

where the prime denotes the derivative with respect to the scale factor a and we choose our initial conditions to be setup 
in the matter dominated era with 5 m (ai) = a,i and S' m (ai) = 1, where Zi = 30 is sufficient. The growth factor is then 
given by D(z) = 5 m {z) / 5 m {z = 0). Furthermore, p m (to) is the present day matter density, f2 m (a) = p m (a)/p cr it(a) is 
the matter density in units of the critical density and om is the over-density with mass M today, where 
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and M = 4:irp m (to)R 3 /3. The window function in (5) is that of a spherical top hat with W{x) — 3(sin(x)/x 3 — 
cos(x)/x 2 ). 

The final missing ingredient is the linear power spectrum P(k), where A 2 (fc) = 4irk 3 P(k). We define the transfer 
function T(k,z) via P(k) = Ak n [T(k, z)} 2 , where n is the spectral index. CMBFAST [78] is used to compute the 
shape of the transfer functions in k space, where it is sufficient to approximate the perturbations of a dark energy 
model with the corresponding ACDM model, since clustering on small scales is not affected by the presence of a dark 
energy component. We fix the baryonic component to have Sl^h 2 = 0.019 in agreement with Big Bang Nucleosynthesis 
constraints [79,80]. To normalize the power spectrum, we use as, the over-density in a spherical region of 8/i _1 Mpc, 
which is traditionally used to quantify the amplitude of the power spectrum on small scales. 

For our fiducial model we choose erg = 0.9 as suggested by recent CMB measurements [81]. Note that this value 
is higher than the one inferred from recent X-ray measurements [82-84]. One might think that this discrepancy 
is not very worrying since as is a parameter which we are trying to determine. However, since the number of 
objects found in a given survey is very sensitive to as and the statistics is Poisson distributed, using a much smaller 
value, as = 0.7 say, would substantially weaken our conclusions. Also we assume that the universe is flat and 
H a = 72kmsec _1 Mpc _1 [85], O m (to) = 0.3 and a spectral index n = 1 [86-88]. For the dark energy equation of state 
parameters we choose u> = —0.8 and w\ = 0.3 which is allowed by current data [90,91,86,87] and is compatible with 
an interesting Supergravity related Quintessence model [61,69]. 

In Fig. 1 we show the comoving mass density for a ACDM cosmology (w=-l). The rapid decrease for large cluster 
masses and large redshifts is clearly visible and we expect a low number of clusters in these ranges. This choice of 
comoving number density is similar to the usual Press-Schechter (PS) formalism [92]. However, the PS approach 
shows some significant differences from the observed and simulated comoving number density. The simulations show 
less low and an increased number of high mass clusters (see also Fig. 13) compared to the PS prescription. A much 
better fit is achieved from using PS if one assumes an aspherical collapse [93,94,41]. Also note that the mass function 
has an uncertainty [41,94,95,77,96] and we will come to this problem when we discuss systematic errors in section V. 

In the following we will discuss the dependence on cosmology of the number of objects per square degree above a 
given mass limit Mn m = 2.35 x 1O 14 M . We will see in section III that Mi; m = Mu m (z), but using a fixed limiting 
mass should provide some intuitive information. In Fig. 2 we present the surface density of clusters (top left) and its 
constituents, the comoving number density, A com = , . -^dM, (lower left), the growth factor (top right) and 

the volume (lower right). As a base for comparison we use the fiducial cosmology described above and then change 
subsequently one parameter. As expected the strongest dependence is on as and Q m (to)- However, we see a weak, but 
nonetheless significant dependence on the equation of state parameters Wg and Wi . The dependence on the spectral 
index n is very weak and hence we fix n — 1.0 in the subsequent analysis. We also see that the strongest dependence 
must come from the linear growth factor (upper right panel), since minute changes lead to significant changes in the 
surface density. This is most obvious from the change in the O m (dashed line) , where the volume element (lower right 
panel) would predict a lower number of clusters, which is balanced by the faster growing linear perturbations (top 
right). This is clear, since the dependence of the surface density on the growth factor is exponential (Eqn. 3). 
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FIG. 1. The comoving mass density dn/dM in units of [Mpc 3 Mo] 1 from simulations of the VIRGO consortium [77] using 
a ACDM cosmology. The lines on the base of the plot represent contours of constant dn/dM. 

III. THE LIMITING MASS 

In the previous section we assumed that the mass limit of the survey is constant. However, in a realistic observing 
situation the survey is limited by the flux and not mass. Hence, we need a relation between the flux limit Su m and 
the mass limit Mn m . We should note that, contrary to what is often explicitly stated or implied in the literature, 
the surveys are not mass limited in the sense that there is single limiting mass which applies across a wide range of 
redshifts. If one were to use a fixed limiting mass, one would be either forced to throw out many objects in order to 
make the survey complete to some high threshold, or put up with lack of completeness at both low and high redshifts 
with a smaller threshold. Here, we will discuss how to compute the limiting mass by first making the extreme, and 
generally incorrect assumption, that the cluster is point-like in the telescope beam. We will then extend our results 
to clusters which are resolved. A more detailed discussion of this problem can be found in the relevant section of the 
companion paper [98] 

A. Principal dependence on cosmology 

We first construct the relation between the virial mass and flux limit of the survey assuming that the source is 
point-like in the telescope beam. The total flux density and the brightness temperature are related by [1,2,97] 

S„ = 2u 2 ATf(x) = 2is 2 T CMB f{x)Y , (6) 

with f(x) — x 2 e x g{x) / (e x — l) 2 and g{x) — x/tanh(izi/2) — 4, where x — 2ttv/Tcmb is the dimensionless frequency. 
The Y parameter is given by the integrated y-distortion Y — J dily(9), with 

y{6) = ^ fdl n e T c = ^t c (9) , (7) 

TO C J TO 

the line of sight integral over the gas pressure in the cluster. The Thomson scattering cross section is or and to c is 
the electron mass. We have introduced the electron density weighted average temperature (T c ) n = J dln c T c / J dln c , 
and r e {9) = <tt / dln e the optical depth. 
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FIG. 2. In the top left we show the surface density of clusters above a mass limit of M lim = 2.35 x 10 14 M Q as a function 
of redshift. The top right is the normalized growth factor D{z), the lower left the comoving number of clusters N coln (z) and 
the lower right is the volume element dV / (dzdQ) . The solid line is the fiducial cosmology (see text), the dotted line a ACDM 
cosmology, the short dashed line is for Q m (to) = 0.4, the long dashed line for n = 0.9 (almost underneath the solid line), the 
dot - short dash line for as = 0.72 and the dot - long dash line for Ho = 65 kmsec _1 Mpc _1 . 



It is conventional to normalize the optical depth in terms of the virial mass of the cluster using 

Mvir/gas = / d 3 r n e (r) = d\ [ dQdln c , (8) 



fi e m p 



from which one can deduce that 



r e (0)=a T ^^^( [MSI) , (9) 
d% fJ, c m p \J J 

where t q (8) = t c (O)£(0), / gas is the intra cluster gas fraction, which we choose / gas = 0.12 throughout this paper [99], 
fi e = 1.143 the mean mass per electron and m p the proton mass. In the general the profile £(0) will depend 
on the virial radius i? V ir and redshift z via the angular diameter distance, which for a flat universe is given by 
<Ja(z) = (1 + z)^ 1 f Q H^ 1 (z') dz' . Using this normalization, it is easy to see from this that the flux density is given 
by 

to ii e m p d\ 

In order to solve (10) for the mass for a given flux limit, we need to know the relation between mass M v j r and 
the cluster temperature (T c ) n . We assume that clusters are virialized objects which are in thermal equilibrium. This 
assumption seems to be confirmed by recent X-ray observations [84] , though it is expected that due to ongoing mergers 
and heat input, this assumption does not hold for high redshifts [100]; an issue we will discuss in section V. If we 
assume that the cluster is virialized the kinetic energy at virialization -E^?" is given by 

E^ = -\u G + U Q , (11) 
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with Ug — — 3GAfy ir /(5i? v ir) the potential energy due to the gravity of the matter. By integrating 6q = pq + 3pq 
over the sphere, with equation of state factor w and assuming that pq is smooth on cluster scales, we obtain the 
general expression for the potential energy 

Uq = m^TrGpq (a)R 2 vir M vir . (12) 

Note that (12) is the generalization of the well known result for a cosmological constant [101], and is different to the 
form used in [40]. We can relate the virial mass of the cluster to the virial radius, by assuming spherical symmetry 
Af vir = 47ri?3. r/9clustcr /3 j where Pciustcr is the mean density of the cluster. The mean kinetic energy E^^} — M v - lr (V^ T ') /2 
is given in terms of the root mean square of the velocity dispersion of the cluster. If the cluster is in thermal equilibrium, 
the virialization temperature is given by [102-104] 

(T e )n = \nm v {VX) , (13) 



with p = 0.59 the mean mass per particle. Hence, we deduce that 



M vir _f(T c ) n \ 3/ \ A ^-1/2 



(A c E(z) 2 Y 



l + (l + 3«,)°# 



-3/2 



(14) 



where T* is a normalization factor which can be deduced to be T* w 0.5 under the assumption that the process of 
virialization only involves gravitational heating . Both simulations [103,41] and observations appear to be somewhat 
inconsistent with this value. For our subsequent analysis we choose T* = 1.6 and we will come back to this in 
section VA. Note that we have expressed the mean cluster density in terms of the over-density A c of the cluster at 
virialization, with Cluster = A cj o crit - 

We can solve numerically for the over-density at virialization by applying the spherical collapse model 
[105,106,38,40]. We follow essentially the discussion in [40] where the background cosmology is described by the 
Friedman equation 

a\ 2 8ttG , . ,„ 

-J = — (fim + pq) ■ (15) 

and the collapsing cluster of radius R by 

R AttG r/ „ „ . . „, 

= — [(1 + 3w) Pq + ^cluster] ■ (16) 

As before, we make the standard assumption that the dark energy component does not clump on the relevant 
scales. We can solve this coupled system numerically and obtain the density of the cluster at turnaround t = t ta , 
Pciustor = CtaPm, an d obtain a best fit for the over-density 

Cta(Zta) « (jfj 2 f7 m (z ta )- - 79 + - 26 ^( z -)- - 08 (— 0-0.21^(1+^) . 

We can then scale this result to the time of virialization, t = t v = 2i ta , with 

A C (Z V ) = Cta(^ta)O m ^) (^) (\^r) 3 , (18) 



where we have introduced 

-Rvir 1 - »7v/2 



i? ta 2 + rh - 3/2r/ v 

the ratio of the radius at virialization (or collapse) of the cluster to the radius at turn around with 



(19) 



Vt = - [1 + 3 W (z ta )] -r%%^ , (20) 



Cta^m(Zv) \1 + Zta. 
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Note that again these expressions differ from [40] because of the generalized potential for the dark energy component 
in (12). We have now all the necessary ingredients to calculate the virial mass of a cluster for a given Sunyaev- 
Zel'dovich flux decrement. However in order to calculate the limiting mass in (3) we need to transform the virial 
mass from an over-density of A c to an over density of A = 200. In order to do this we assume that the matter in 
the cluster is distributed according to a Navarro- Frenk- White (NFW) profile [107], with a concentration parameter 
c = 5. This allows us to rescale the cluster mass from M v ; r to M2oo- In Fig. 3 we illustrate the limiting mass and 
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FIG. 3. The limiting mass (left) and surface density of clusters (right) in the point source approximation for a survey with 
limiting flux density S\i m = 4.4 mjy and frequency v — 150 GHz. The different lines correspond to different cosmological 
parameters with the key the same as in Fig. 2. 



redshift distribution of clusters using the point source approximation for the same cosmologies as in Fig. 2. First, 
notice that the limiting mass depends on the cosmological parameters. As well as the standard dependence on Hq, we 
see dependence on J7 m and wq which would be ignored by assuming that the limiting mass is constant. Note, that the 
dependence of Mn m on wq now lowers the overall number of clusters (dotted line) , balancing the slight increase of the 
number of clusters for the ACDM universe, compared to the fiducial model in the case of a fixed limiting mass (see 
Fig. 2) . For an increasing Q m the limiting mass is decreasing resulting in an even more enhanced number of clusters 
compared to the fiducial model. Moreover, in this approximation the limiting mass goes to zero at low redshifts. 
Although this will be seen to be unphysical when telescope beams are small, it will be a good approximation when 
the beam is large. The dependence of the limiting mass on the cosmological parameters feeds into the predictions for 
dN/dz, but the effects of parameters such as as and n is unchanged from that for a constant limiting mass. 



B. Extended clusters and survey beam 

So far we have not accounted for the possibility of clusters being extended objects within the telescope beam which 
could lead to a loss of flux. If there is a beam B(9) present, integrated y-distortion in (6) becomes [108] 

where we assume a Gaussian profile of the beam, with B(6) = exp[— 8 2 / (2of )] and 07, is related to the full-width-half- 
maximum (FWHM) of the beam by db = 0FWHM/V8 In 2. The equation for the observed flux for a beam directly 
over the centre of the cluster is then [98] 

^ = 2^T CMB (T c ) n f gas M vil 1 c) (23) 

m c ^ c m p d £ K 

with 

jdnB(e)qe) 

[ j<mB(e)jdnc(e)- [ ' 
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In general this function will depend on M v j r via i? V i r , so a simple analytical solution for M v i r in terms of Sn m and 
z is not possible and therefore we must solve (23) numerically using the Newton-Raphson method. The expression 
(23) is the flux in a single beam, however, if the cluster is larger than one beam area , beams can be combined with 
a resultant increase in the noise; a process often known as beam smoothing. The optimum number of beams is given 
by 

as discussed in [98]. In reality iVb must be an integer, but we will allow it here to take any value greater than one. It 
was shown in [98] that this does not unduly exaggerate the predicted number of clusters detected. 

In most situations there will be a limit to the extent to which one can beam smooth, for example, due to the 
primary CMB anisotropy. In the subsequent discussion we will assume that 0fwhm = 20'. Note that for multi- 
frequency observations like Planck it may be possible to effectively filter out the SZ signal due its distinct spectral 
signature. This in principle would allow to combine even more pixels, however we do not include this effect in our 
analysis and we do not expect it to be significant, since it will be important only for objects at low redshifts. 

We assume that the cluster is in hydrostatic equilibrium and can be described by a truncated isothermal (3- 
model [109,110], with an electron density distribution of 

n c (r) = \ n ° i 1 + lj ) 2 r< R ™ f (26) 
r > i? vir 

where R c is the core radius of the cluster defined in terms of the virial radius R c = R vn /a. The projected profile 
function is then 



a2\ i-iP J 

t(0)= (i + 



( a2 - d y e . 2 A 1/2 r 

\ 1+5*755 ) >P 



J{a,(3) 



(27) 



where J(a, b) = f (l + x 2 ) 3b ^ 2 dx and 9 C = R^/dx- For the purposes of our discussion we will assume that (3 — 2/3 
and a = 10 [108], in which case 

J(a, 2/3) = tan" 1 (a). (28) 
Fig. 4 illustrates the limiting mass and redshift distribution of clusters if we use a beam of #fwhm = 1-3'. The 




FIG. 4. The limiting mass (left) and surface density of clusters (right) for a survey with limiting flux Su m = 2.62 mjy per 
beam, a frequency v — 150 GHz and a beam width of #fwhm = 1-3' for an optimal combination of pixels up to a angular size 
of 0pwHM = 20'. The different lines correspond to different cosmological parameters with the key the same as in Fig. 2. The 
light solid line corresponds to the no beam smoothing case, when we fix iVb = 1. 
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cosmology dependence is very similar to the point source approximation. We also include the result if we allow no 
beam smoothing (light, solid line). We see that below a redshift of z = 1 the possibility of beam smoothing becomes 
an important feature. With #fwhm = 1-3' and $fwhm = 20' it is possible to combine up to 60 beams, which leads 
for z = 0.4 to a limiting mass which is 5 x 1O 13 M lower than without smoothing. Furthermore we note that the 
limiting mass is not dropping to zero as it is the case in the point source approximation. For the surface density we 
see how combining the pixels improves the situation in the low redshift region. We have now all the ingredients to 
perform a likelihood analysis for different types of the surveys, which we will discuss in the next section. 

IV. LIKELIHOOD ANALYSIS OF PROPOSED SURVEYS 
A. Proposed surveys 

In order to be able to make some firm statements as to the extent to which one might be able to constrain 
cosmological parameters we have to specify the flux limit Su m , the angular coverage Af2, the frequency of observation 
v and the beam size #fwhm of the survey. This is a large parameter space and it would be prohibitive to explore it 
fully. Hence, we have decided to devote our attention to modelling six instruments which have been proposed and 
which cover a wide range of possibilities. They are the VSA [29], BOLOCAM [16,17], the Arcminute Microkelvin 
Imager (AMI) [18,19], the One-Centimetre Receiver Array (OCRA) [24], the Planck Surveyor Satellite and the South 
Pole Telescope (SPT) [25,26,111]. We have not specifically include the Sunyaev-Zel'dovich Array (SZA) [20-22] and 
the Array for Microwave Background Anisotropy (AMIBA) [23] in our analysis since they are likely to of very similar 
sensitivity to AMI. Other experiments are proposed, but we do not have sufficient information to include them in our 
calculations. We should caution the reader that these instruments use a wide range of different techniques, have/will 
cost wildly different amounts of money and will become operational over a wide range of timescales. It is not our 
intention to discuss the merits of particular experiments rather give a picture of the likely progression of the subject. 

The values which we have chosen to represent each instrument are presented in Table I. The particular values for 
Siim and Ail use a quoted instantaneous sensitivity (either for the whole instrument or individual beams) and the 
optimal depth of the survey using the procedure discussed in [98] . The angular coverage can then be computed from 
the instantaneous field of view and assuming one full year integration time. Very different results would be achieved 
by trading off sensitivity for area according to Sii m /V AO. However, any changes from the quoted values would lead 
to fewer objects being found and would, therefore, be likely to lead to weaker constraints on cosmological parameters. 
This statement is only true for our particular modelling of the distribution of clusters. 

It is worth making some comments on how we have done this for each of the different instruments: 

VSA : A 14-element interferometer based on Tenerife. We use the specification for a proposed upgrade which should 
have an instantaneous sensitivity of / = 0.5 Jy sec 1 / 2 on a field of view of ^fov = 1 deg 2 . 

BOLOCAM : This instrument has 144 instantaneous beams and is being used on the 8m telescope at the Caltech 
Sub-millimetre Observatory on Hawaii. We estimate #fwhm = 0.8' at 150 GHz and each of the beams has / = 
35 mJy sec 1 / 2 . Atmospheric emission could well require some beam switching, but for our calculations we assume that 
all 144 beams measure the sky providing an instantaneous field of view Ofov = 0.025 deg 2 . 

AMI : Consists of the eight 13m dishes of the Ryle Telescope and a new array of ten 3.7m antennas. It should 
have V^fov = 21' and / = 20mJysec 1 ' 2 . The extra resolution provided by the Ryle baselines should allow for an 
investigation of the potential problem of radio point sources inside the clusters which would dilute the signal at low 
frequencies. 

OCRA : Currently funded is a 10 beam system which will use 30 GHz receivers similar to those developed for the 
Planck LFI system. This will be mounted on the 32m telescope at Torun in Poland. Since the atmosphere is likely to 
be significant at this location the receivers will be used in pairs each with / = SmJysec 1 / 2 . We present results for a 
proposed 100 beam system which should have ^fov = 0.034deg 2 . 

PLANCK Surveyor : Here, we use the likely sensitivity for 2 whole scans of the sky which should take 14 months for 
the 100 GHz channel on the HFI instrument. We conservatively assume that it should be possible to extract the SZ 
signal over 1/2 of the sky. We have not attempted to optimize the yield of this survey since the satellite will cover 
the whole sky. 

SPT : This is an ambitious project to build a 10m telescope at the South Pole and mount on it a 1000 element 
bolometric array. It should be possible for such an array to achieve Sum = 5mJy on A£l = 4000 deg 2 [111] 
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The number of clusters we predict for the fiducial cosmology and the setups of Table I are w 60 for VSA, ss 400 for 
BOLOCAM, « 170 for AMI, w 460 for OCRA, w 4500 for SPT and « 5200 for the Planck surveyor. 





Sum [mjy] 


v [GHz] 


^fwhm 


AO. [deg 2 ] 


VSA 


5.75 


30 


8.0' 


300 


BOLOCAM 


1.64 


150 


0.8' 


130 


AMI 


0.58 


15 


4.5' 


230 


OCRA 


0.30 


30 


1.1' 


140 


Planck 


41.5 


100 


9.2' 


20600 


SPT 


2.62 


150 


1.3' 


1430 



TABLE I. The experimental parameters for the six different SZ surveys assuming approximately one year integration time. 




0.5 1 1.5 2 0.5 1 1.5 2 

z z 

FIG. 5. Limiting mass (left) and number of clusters (right) in redshift bins of Az — 0.1 for the different experimental setups 
from Table I. The dot - short dash line is for BOLOCAM, the dot - long dash line for VSA, the solid line is for the AMI survey, 
the long dashed line for OCRA, short dashed lines for the 100, 143 and 353 GHz channels of the Planck Surveyor (note that 
we only use the 100 GHz channel in the further analysis since it predicts the largest number of clusters) and the dotted line for 
SPT. The data points are from a randomly generated sample assuming Poisson errors. Note that we have approximated the 
errorbars in the plot with a Gaussian distribution of width y/ AN. The plots are all for the fiducial cosmology. 



In Fig. 5, we present the computed mass limit of the proposed surveys. Due to the fact that we have optimized the 
survey yield, instruments with similar beam sizes have similar limiting masses, that is SPT, OCRA and BOLOCAM 
have the lowest limiting masses and are approximately the same, as are VSA and Planck which have the highest 
limiting mass because of their large beams and low sensitivity per beam. On the right in Fig. 5 we show the surface 
density of clusters in bins of Az = 0.1. VSA should find a handful of clusters in each bin while AMI could discover 
up to ~ 20 clusters in some bins around the peak near z = 0.3. BOLOCAM and OCRA will observe a few tens of 
clusters in the redshift bins around the peak region and many more at high redshifts than Planck. This is is due to 
the small beam size. Although SPT will observe fewer clusters than Planck, in the region above redshift z = 0.5 the 
cluster yield of SPT is far larger than the one from Planck, again due to its comparably small beam. 



B. Poisson statistics 



In the following we will discuss the ability of these experiments to constrain cosmological parameters. Therefore 
we need to define a statistic to treat this problem quantitatively. We assume that the measurement is dominated by 
Poisson statistics. If we measure Nt,in,i clusters in a particular redshift bin Zi the overall probability to observe {fii} 
clusters in JVdat bins is then 

JVdat J\J n i 

P ({ni} | { A bin .J) = [J __2ii e -"w».< . (29) 
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We would like to test how probable it is that a given theory fits our measurement {Abi n ,i, Zi}. Therefore, we have 
to calculate p \\N^ aTy {z%)\ | {iVbin.i}) , with N^ OTy (zi) = /J!'^/2 dN/dz(z) dz. We can build a log-likelihood 
statistic by defining 

C = -2 lnp = ~2J2 <„° 0ry (^) IniVbm^ - N hirl<i - \nN^( Zi ) , (30) 
i=i 

which is known in the literature as Cash C statistic [112]. We will sample the log-likelihood function in (30) around 
our fiducial cosmology, where we assume that we can get redshifts out to z max = 1.5 with an accuracy of Az — 0.1. 
This should be feasible with photometric redshifts from instruments like the SDSS and VISTA [74-76] . For an efficient 
way of sampling we adopt a Markov chain Monte Carlo sampling method [113,114,87], where we typically compute 
half a million of samples to achieve convergence. For the likelihood scanning we vary Hq, f2 m , erg, wo and Wi and 
keep the other parameters fixed. In order to speed up the likelihood calculation we also approximate the calculation 
of cr(R) in (5) by the expressions given in [115,116]. 




C. Mock likelihood analysis 



In Fig. 6 we present the joint 1 — a likelihood contours for the mock surveys. We have marginalized over all the 
other parameters and the only prior assumption made was that Hq — 72 ± 8 kmscc _1 Mpc _1 as suggested by the 
Hubble Space Telescope Key Project [85]. Note that we include a larger number of parameters in our analysis than 
in previous analyses which naturally tends to increase the size of the errorbars. It is clear that the results from VSA, 
BOLOCAM, AMI and OCRA are unlikely to improve on current constraints on the cosmological parameters, but 
that they will provide extra independent, complementary information which can be used in conjunction with, for 
example, that from the primary CMB. It is possible that these surveys could be used in conjunction with X-ray or 
weak lensing observations of clusters to probe the gas physics of clusters and as well as to constrain cosmological 
parameters, providing further insight for future more powerful surveys. 

Our analysis shows that Planck and SPT are likely to yield tight constraints on erg and Q m . They will be also able 
to start to constrain the equation of state parameters. However, it seems that the surface density of clusters alone will 
not be able to distinguish the fiducial model from a ACDM cosmology. In order to achieve tighter constraints on the 
equation of state it will necessary to use additional complementary information from SNe observations [117], baryon 
mass fraction measurements [88] or other cluster abundance measurements [84] . In the subsequent discussion we will 
focus on the SPT survey. 

In Fig. 7, we investigate the use of prior information on f2 m and erg. First, we include an uncertainty of Af2 m = ±0.04 
(solid). This accuracy has been already achieved for measurements of the baryonic gas fraction in clusters using X-ray 
techniques [88]. It was established that this uncertainty is relatively insensitive to the inclusion of a dark energy 
component and hence we use it without modification [89]. We see that such a prior does not significantly improve 
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FIG. 7. The 1 — a joint likelihood in the W\ — u>o plane for the SPT survey. On the left, standard SPT (dotted), with a prior 
on il m = 0.3 ± 0.04 (dashed filled contour) and with prior on Q m = 0.3 ± 0.02 and as ~ 0.9 ± 0.01 (inner dot-dashed filled 
contour). The dashed line is for the Supernovae Acceleration Probe - SNAP. On the right, z max = 1.0 (solid), standard SPT 
(dotted), SPT with Az = 0.01 (dot-dashed). 

our ability to constrain the dark energy parameters. In order to improve on this, we have also put very tight priors 
on f2 m and as, with AS! m = 0.02 and Aug = ±0.01 (dot-dashed); a level of constraint which could be possible using 
X-ray and lensing observations [84] . This serves illustrate how sensitive the cluster abundance measurement is to the 
value of as since the constraint on the dark energy is substantially improved. If there is a tight prior on a% then other 
cosmological parameters can be tightly constrained using SZ surveys. 

The dashed open contour is the 1 — a likelihood for the proposed SNAP survey. We assumed for this analysis 
that SNAP will discover about 2000 SNe out to redshift of z = 1.7 [68,69]. We see that SZ surveys provide valuable 
complementary information to SNe surveys. The SNe seem to constrain the constant part of the equation of state, 
wo, much more tightly than the evolving part w\, while SZ clusters appear to constrain both parameters equally. 
This is because the surface density of clusters essentially constrains cosmological parameters via the linear growth 
factor D(z) and the volume element dV/(dzdd) t while SNe constrains on the magnitude-redshift relation. As a note 
of caution we find that even using SPT with the tightest priors as well as SNAP one cannot establish the evolution 
of the equation of state of the fiducial model. However, if we marginalize over w\ , a difference between cosmological 
constant and wq > — 1 can be established albeit with limited significance. 

In Fig. 7 (right) we illustrate how our knowledge of the cluster redshifts can affect our conclusions. The outermost 
solid line likelihood contour is for the SPT survey with z max = 1. In this case, the errorbars are much larger than 
those expected when we assume the measurement of redshifts out to z max = 1.5. The dashed contour is obtained if 
we have redshift bins of Az — 0.01 instead of Az = 0.1. We see that there is very little difference between the error 
contours for these two different binning values suggesting that photometric redshifts should be sufficient at least for 
this application. This can be seen by looking back at Fig. 5 where we see that Az — 0.1 seems sufficient to map the 
shape of the surface density of clusters since it evolves much more slowly with redshift. 

V. STATISTICAL AND SYSTEMATIC UNCERTAINTIES 

A. Mass - temperature relation 

In order to calculate the surface density of clusters above a certain mass limit it is necessary assume a mass- 
temperature relation (14). The derivation of the mass temperature relation assumed that clusters are completely 
virialized objects in thermal equilibrium governed only by gravitational physics. However on-going mergers, incomplete 
virialization and numerous possible non-gravitational heating mechanisms, such as AGNS and SNe, can modify this 
assumption. We will, therefore, discuss the consequences of modifications to the mass-temperature relation. First 
we will concentrate on the overall normalization amplitude T„ which we have already pointed out is deduced by 
normalization to the results of numerical simulations. From the compilation of different simulations and observations 
in [123] we see that T» can be anywhere in the range T* = 1.0 — 2.2. We have chosen T* = 1.6 for our analysis which 
is in the middle of this region, although not preferred by a particular measurement or simulation. 
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FIG. 8. The surface density for the SPT survey for different values of T*. The dotted line is for T* = 1.0, the long dashed line 
(partially hidden) for T» = 1.44 (-10%), the dot-dashed line for T» = 1.76 (+10%) and the dashed line for T, = 2.0. The middle 
solid line is for the fiducial model, the upper solid line for f2 m = 0.32 and the lower solid line for Ho = 65 kmsec -1 Mpc -1 . 
Clearly, differences of ±10% are of the same order of magnitude as the constraints on cosmological parameters which we are 
trying to establish. 



In Fig. 8, we show the variation of the redshift distribution of clusters with varying T*. We also show the values for 
the fiducial model and models with f2 m = 0.32 and Ho = 65 km sec Mpc -1 . We see from this that an uncertainty in T* 
changes the amplitude in the surface density but not the shape and, therefore, we expect that the degeneracy between 
the equation of state parameters, which mainly defines the shape of the curves (see Fig. 4), is largely unaffected by 
this uncertainty. However, we expect degeneracies with as and Q m . 

In order to test this, we first study the bias we introduce if we analyse mock data, created using a different value 
of T» than the one we used to perform the analysis. In the top two and lower left panel of Fig. 9 we show the results 
of this. The dotted line corresponds to a model created with T* = 1.86 and the dashed line to T* = 1.42. In all cases 
we performed the analysis for the fiducial value T» = 1.6. This roughly corresponds to a bias in the limiting mass 
of ±20%. We see clearly the bias introduced in the as and f2 m measurements by choosing the wrong value of T*. 
However, the effect on the dark energy parameters wq and w\ (Fig. 9, lower left) seems to introduce only a slight extra 
bias toward larger wq values and only increases or lowers the errorbars corresponding to the change in the overall 
number of observed clusters. 

This becomes even clearer if we include the parameter in the likelihood analysis and marginalize over T*. The 
result of including T» is presented in Fig. 9, lower right panel. We see that this only increases the errorbars in the wq 
direction. This is because wi (and also partly wq) influence mainly the shape of the surface density and not so much 
its amplitude, as suggested above. 

So far we have assumed that the mass-temperature relation applies universally. However, at this stage it is only 
an observed correlation and is likely to have considerable scatter. Next, we investigate the inclusion of statistical 
uncertainties in the limiting mass. We incorporate this via a "selection function" <j>(M,z) [25], 



o 

with 



oo 



dN . „ dV 

dz 



M — Mii m (z) 

SM lim (z) 



(32) 



where 5 is the relative statistical uncertainty in the mass limit Mu m . In Fig. 10 we illustrate the influence of different 
values of S on the observed surface density. A general result is that the inclusion of this effect using a symmetric 
selection function will always increase the overall number of observed clusters since there are many more objects just 
below the mass threshold than just above it. We see from the figure that an unknown uncertainty in the limiting 
mass will also increase the uncertainty in constraints on the cosmological parameters, but if it is known and well 
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FIG. 9. The 1 — a likelihood contours for the SPT survey. In the top two panels and the lower left panel we show the effect 
of the mass-temperature normalization bias, if we assume T* is 16% larger (dotted) and 12% lower (dashed) than the fiducial 
value of T, = 1.6. In the lower right panel T» is included as a free parameter. The inner light shaded contour is with fixed 
T* = l.G and the outer contour marginalized over T, . 
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FIG. 10. The surface density of clusters for the SPT survey. The lowest solid line is for the fiducial model, the middle solid 
line for f2 m = 0.32 and the upper solid line for as = 0.95. The dotted line is for a statistical uncertainty in the limiting mass of 
S — 20%, the short dashed line for 5 — 30% and the long dashed line for 8 — 40%. The data points are from a simulated mock 
catalogue for the SPT survey. 
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understood it should be possible to incorporate it into the analysis. Since the changes in the surface density due to 
S only become larger than the expected Poisson errors for SPT when the scatter is greater than 6 — 15%, we 
estimate that a statistical uncertainty below this level should be acceptable for this survey. The Poisson errors for 
the other surveys are much larger and hence it should be possible for them to withstand much larger values of S and 
still give sensible constraints. 

The final uncertainties in the mass temperature relation which we shall investigate are a different power law depen- 
dence between M and T, and a different overall dependence on redshift [100,25]. We re-define the mass temperature 
relation to be 



(T e Y = T, (A c E(zf 



,1/3 



1 + (1 + 3w 



n Q (z) 



(i 



M v 



(33) 



where if e = 1 and £ = 3/2 we obtain the standard mass temperature relation (14). Note that new values for e and 
£ would require a recalibration of the relation and hence a different value of T*. However, since we will only restrict 
ourselves to a qualitative discussion we will not incorporate this effect. The parameter £ could model non-gravitational 
heat input. Since smaller clusters (that is, groups) are preferentially affected by these processes, we expect £ > 1.5. 
Observations and simulations suggest values of 1.48 < £ < 1.98 [118,99,119-121]. The parameter e models deviations 
from complete virialization. On one hand, clusters at early times might not be completely virialized, hence e < 1. 
However, clusters which have ongoing mergers or some other form of energy injection could be much hotter than 
expected and hence have e > 1. At this stage there is no observational preference for any particular value of e. 
In Fig. 11 we illustrate the surface density for the SPT survey for different values of £ and e. We allow £ to vary 
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FIG. 11. The surface density of clusters for the SPT survey with changes in the mass - temperature relation. The lowest 
solid line is for the fiducial model, the middle solid line for f2 m = 0.32 and the upper solid line for as = 0.95. The dotted lines 
are for a change in the power law between limiting mass and temperature. The lowest dotted line is for £ = 1.6, the middle 
dotted line for £ = 1.7 and the top for £ = 1.9. The long dashed lines are for changes in the redshift dependence of the M — T 
relation. Form top to bottom the long dashed lines correspond to e = 1.5, 1.2, 0.8, 0.5. The data points are from a simulated 
mock catalogue for the SPT survey. 



in the range £ = 1.6, 1.7, 1.9 and see that we observe more clusters for higher values of £. This is clear since from 
(10) we obtain M vlr cx S$ /{£+1 \ This again will increase the uncertainties mainly in the parameters which define 
the amplitude of the redshift distribution, like O m or erg. We vary e between 0.5 and 1.5. A recent Fisher matrix 
analysis of SZ cluster counts for the SPT survey obtained uncertainties of A£ w 0.0064 and Ae « 0.46 [25], if they 
are included as parameters. Although we believe that a Fisher matrix analysis can only give crude errorbars, we 
expect degeneracies to not restrict our ability to extract the cosmological parameters. We should note that follow 
up measurements which constrain the mass of the cluster independently using, for example, X-ray measurements or 
weak lensing, could constrain the non-standard mass temperature relation considerably [25,100] and, hence, improve 
the veracity of constraints on the cosmological parameters deduce from SZ surveys. 
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B. Mass - temperature relation and as 



Recent years have seen a wide range of values reported in the literature for the power spectrum normalization 
cr 8 . Observations with weak lensing, X-rays, Large-Scale-Structure and the large scale SZ effect suggest values in the 
range 0.61 < cr 8 < 1.05 [122]. If interpreting a particular observations requires the conversion from temperature to 
mass, as is the case for the SZ effect and X-ray observations, then it is important to know the value of T» which has 
been used, since there is a direct degeneracy between cx 8 and the normalization with f2^ 6 <jg oc T~ ' 8 [122]. As 
noted in the introduction one problem of our analysis is that the choice of a lower value for as in our fiducial model 
would have resulted in a lower number of observable clusters and hence much weaker constraints on the cosmological 
parameters due to the Poisson statistics. However, the analyses of X-ray measurements which result in 0.7 < erg < 0.8 
use a larger value of T* [96,71,122]. We will now discuss various aspects of this issue. In Fig. 12 we illustrate the 
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FIG. 12. The surface density of clusters for the SPT survey. The solid line is for the fiducial parameters and set up. The 
dotted line is for the same setup, but as = 0.72. The dashed lines are for the as = 0.72 with T* increasing from 1.8 (lowest 
dashed line) via T* = 2.5 (middle dashed line) to T„ = 3.5 (top dashed line). 

consequences of combining lower values of as with larger values of T*. The solid line is for the fiducial model and 
the SPT survey. The dotted line is for the same set up with as — 0.72, which seems to be in the preferred range of 
X-ray and Large-Scale-Structure observations [87,88,84]. We sec that the overall number of clusters is considerably 
less for the lower value of as- However, X-ray observations seem to require a lower calibration of the mass temperature 
relation [84,122], which corresponds to an increased value of T*. In Fig. 12 we show the increasing number of clusters 
observed for as T* increases. Still for the marginally maximal acceptable value of T* = 2.5 [123] the number of clusters 
is about a factor of 1.7 below the number observed in a universe with ag — 0.9. Nonetheless it is clear that the lower 
values of as suggested by some lead to the very small number of clusters one would naively assume. 

In order to obtain as many clusters predicted in the fiducial model, but using as = 0.72 it is necessary to increase 
to = 3.5, which corresponds to a 75% lower normalization of the mass - temperature relation, outside the range 
suggested by simulations and measured in reality [123]. Hence, if one is drawn to conclusions of these particular 
observations it seems likely that our fiducial model over-predicts the number of clusters by at least a factor of two 
and the errorbars on the parameters we have predict are optimistic. However, there are observations which suggest 
the opposite. Preliminary observations of secondary CMB anisotropics of the large-scale SZ effect using the CBI 
instrument [124] indicate a value of as > 1 [125]. Although this results remains to be confirmed, it indicates that the 
issue of the value of as is far from settled. We should note that it might be possible to measure the normalization of the 
mass - temperature relation directly with the combination of cluster abundances and weak lensing observations [123] 
and constrain both as and directly. Furthermore, a measurement of the flux of each cluster will enable one 
to perform an internal calibration of the sample and alleviate the uncertainty in the normalization of the mass - 
temperature relation [126]. 
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C. Mass function 



A further uncertainty we need to consider is the mass function itself. In our analysis we used the 2002 results from 
the VIRGO consortium [77] . In order to obtain insight into the uncertainty of the mass function we have compared 
our results to the previously released mass function from 2001 [94] 

*L (Mfl . -0.3 1 6^^ J Lexp { -|0.67-.og [D(2M |-} , (34) 

with the mass defined by an overdensity of A = 324 relative to the matter density which corresponds to an overdensity 
of A = 97 relative to the critical density for the fiducial universe. We have also compared our results with the standard 
PS mass function [92] 

dn , » n [2p m (t ) S c da M ( 5 c \ / QK -> 

dM (Z ' M) = ~V* M ^(zK, dM eXP {- 2D{z)*o*J (35) 
with 5 C = 1.686, where we ignore the weak cosmology dependence of S c in the present discussion. In Fig. 13, we 




FIG. 13. On the left the different mass functions at z — 0. The dashed line is for the Jenkins 2001 mass function [94], the 
solid line for the Evrard 2002 mass function [77] both from the VIRGO consortium. The dotted line is that expected from the 
PS formalism for a mass definition of M2oo- On the right the expected surface density of clusters for these mass functions for 
the SPT survey. 

show the different mass function at z = for the fiducial cosmology. The solid line is for the 2002 mass function 
of the VIRGO consortium [77], the dashed line for 2001 [94] and the dotted line for the PS mass function. We will 
investigate the consequences of these three mass functions for expected yield of the SPT survey. Since SPT has a limit 
of roughly 1.7 x 1O 14 /i -1 M0 we have to consider the mass function above this limit. The PS function dominates over 
both the VIRGO 2001 and VIRGO 2002 mass function up to masses of 10 16 /i _1 M© and then PS mass function drops 
below both VIRGO mass functions. The VIRGO 2001 function is marginally above the VIRGO 2002 function in the 
range below 10 16 h~ 1 MQ. However, above this region the number of clusters is so sparse that the contribution to the 
surface density of clusters is not significant. We also present the surface density for the SPT survey in Fig. 13. We 
see that the PS mass function results in about twice as many clusters and the VIRGO 2001 about 25% more clusters 
than the VIRGO 2002 mass function. It is clear from this analysis, that a precise convergence of the mass function 
is required to constrain cosmological parameters using the surface density of clusters, providing further impetus for 
an extended program of numerical simulations. Clearly the uncertainty in the mass function is degenerate with other 
uncertainties in the mass - temperature relation and o%. However, combining SZ, X-ray and weak lensing of clusters 
could help constrain the mass function as well. 



VI. DISCUSSION AND CONCLUSIONS 



In this paper we have introduced a realistic model for the surface density of clusters that would be observed by 
an SZ instrument including effects of the beam of the survey and the profile of the cluster. In the Section IV we 
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analysed constraints on cosmological parameters that might be expected from such surveys. We found that SPT 
and Planck will be particularly powerful in constraining the standard parameters as and f2 m , and can be used in 
conjunction with other observations, for example SNe surveys, to constrain the dark energy. These surveys can be 
expected on a timescale of around 5 to 7 years. In the meantime less powerful surveys such as those possible using 
VSA, BOLOCAM, AMI and OCRA will provide useful complementary constraints, while allowing the study of the 
clusters themselves, in particular using X-ray follow up. 

We have shown in Fig. 7 that it is sufficient to use redshift bins of Az = 0.1 to constrain the cosmological parameters. 
Essentially the shape of the redshift distribution is mapped out adequately with this bin width. It should be feasible 
for surveys like SDSS or VISTA [74-76] to measure redshifts of this accuracy, although one will have to be careful to 
avoid introducing extra selection effects. 

While Planck, SPT and OCRA will constrain the standard cosmological parameters significantly, it will be difficult 
for them to provide information about the equation of state of the dark energy component. While they do not 
constrain it by themselves, they are complementary to SNe observations as illustrated in Fig. 7. SZ surveys essentially 
constrain f2 m and provide a prior for SNe surveys like the SNAP mission. However, they become considerably more 
powerful if one includes a tight prior on the amplitude of the density fluctuations as, which could be provided from 
complementary large scale structure observations. 




FIG. 14. The probability for wo marginalized over all the other parameters. The solid line is for the OCRA survey, the short 
dot-dashed line for SPT, the dashed line for SPT including the parameter T* in the likelihood analysis and the dotted line for 
SPT with tight priors on fi m and as- 

This leaves us with the question of whether any of the proposed SZ surveys might be able to distinguish the 
fiducial cosmology from a ACDM universe. In Fig. 14 we show the one dimensional likelihoods for the constant part 
of the equation of state wq. We see that the OCRA survey can not constrain the equation of state. However the 
marginalized likelihood of SPT narrows in on the fiducial value of wq = —0.8, while with tight priors on il m and as 
we can distinguish the fiducial model with w = —0.8 + 0.3z from a ACDM cosmology. 

We have also discussed in detail systematic and statistical uncertainties which might affect more ability to affect 
our ability to extract cosmological parameters. The main sources of uncertainty are the mass-temperature relation 
and the mass function; clearly much more work is needed to pin these two quantities down. The inclusion of the 
normalization factor T* in the likelihood analysis only marginally increases the errorbars on the equation of state 
state as shown in Fig. 14. A potentially more significant problem is a modified evolution or power law of the mass 
temperature relation. The Fisher matrix analysis presented in [25] suggests that this problem is not too severe. This 
is because the cosmological dependence of the surface density is mainly due to the growth factor, rather than that 
from the limiting mass. 

The surveys discussed in this paper will provide much more information than just the surface density. In principle 
it should be possible to use the flux information as well as the surface density to constrain cosmological parameters 
as hinted in [100,126]. However, using the flux as well as redshift information is likely to be more susceptible to 
systematic error. This subject is currently under investigation. 

Finally, we should note that it is difficult to strongly constrain the equation of state of the dark energy component, 
and particularly its redshift evolution, using any single method. It is clear that a joint effort will be required with a 
combination of measurements such as those from the cosmic microwave background, large scale structure, SNe, weak 
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lensing and clusters being necessary to achieve this goal. The results presented in this paper show that the SZ effect 
is likely to play a strong role in this over the next decade. 
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